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I Abstract 

We present an unfolding of the codimension-two scenario of the simulta- 
neous occurrence of a discontinuous bifurcation and an Andronov-Hopf bi- 
furcation in a piecewise-smooth, continuous system of autonomous ordinary 
differential equations in the plane. We find the Hopf cycle undergoes a grazing 
] bifurcation that may be very shortly followed by a saddle-node bifurcation of 

the orbit. We derive scaling laws for the bifurcation curves that emanate from 
the codimension-two bifurcation. 

o 

O ■ 1 Introduction 



A system of differential equations is said to be piecewise-smooth if it is everywhere 
smooth except on some codimension-one boundaries called switching manifolds. This 
paper is concerned with such systems that are continuous everywhere but non- 
differentiable on switching manifolds. Piecewise-smooth continuous systems have 
been utilized to model a range of diverse physical situations, for instance vibro- 
impacting mechanical systems [H E], switching in electrical circuits [HI HJ [5] and 
various non-smooth phenomena in biology and physiology P, [7j . 

The interaction of invariant sets with switching manifolds often produces bifur- 
cations not seen in smooth systems. The last two decades have seen an explosion 
of interest in such bifurcations and many new results, see for instance O HI H] and 
references within. In the neighborhood of a single switching manifold a piecewise- 
smooth, continuous system may be written as 

r FW(x), Hix)<0 
1 F^^Hx), H(x) > ' 



*D. J. W. Simpson and J. D. Meiss gratefully acknowledge support from NSF grant DMS- 
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where x E , p(L)^p(R) ; ^ ^jv Qk and : ^ M is sufficiently 
smooth. The switching manifold is the set S = {x E | H{x) = 0} and by 
continuity, F^^) = F^^'> on S. 

A point x* E is an equilibrium of the left-half-system if F^^\x*) = and 
is said to be admissible if H{x*) < and virtual otherwise (and vice- versa for 
the right-half-system). By continuity of ([T]), an equilibrium that lies exactly on S 
is an admissible equilibrium of both smooth half-systems. This co dimension-one 
phenomenon generally gives rise to what is known as a discontinuous bifurcation 
(sometimes called a boundary equilibrium bifurcation) for which there are two basic 
generic scenarios. Either equilibria in the two half-systems coexist, collide and 
annihilate on the switching manifold {non-smooth saddle-node) or one equilibrium is 
admissible on each side of the bifurcation [persistence) . In addition, other invariant 
sets may be created at the bifurcation, such as a periodic orbit in a manner akin to 
an Andronov-Hopf bifurcation [9l [TOj . 

Even though linear terms of an appropriate series expansion of the system gener- 
ally completely determine dynamical behavior local to a discontinuous bifurcation, in 
a general A^-dimensional system the bifurcation may be extremely complex. Much 
work has been done for low [N < 3) dimensional systems Il2l [HE]. In two- 
dimensions, if a periodic orbit emanates from a generic discontinuous bifurcation, 
it must encircle an equilibrium of focus type ^H]. In this paper we unfold about 
the CO dimension-two point where the focus-type equilibrium has purely imaginary 
eigenvalues at the crossing point. This scenario corresponds to the simultaneous oc- 
currence of a discontinuous bifurcation and a (smooth) Andronov-Hopf bifurcation 
and has recently been observed in a model of yeast growth [13]. 

We find that as parameters are changed to move away from the Hopf bifurcation, 
the associated Hopf cycle grows in size as usual, until grazing the switching manifold. 
No bifurcation occurs at the grazing point in the sense that the phase portrait 
does not change topologically (because the system is continuous at the switching 
manifold [S]). However, very shortly beyond the grazing bifurcation a saddle- node 
bifurcation of the orbit may occur. We derive a condition governing when this occurs 
and a scaling law describing how close to the grazing the saddle-node bifurcation 
occurs. Our results are presented in Theorem [H 

The remainder of this paper is organized as follows. In [|2] we transform the 
general system to a normal form involving companion matrices and state the the- 
orem. In ^we provide a simple example and use it to illustrate the theorem and 
numerically verify predicted scaling laws. §1] outlines our method of proof and ^ 
presents a proof of the theorem. Conclusions are presented in ^ 

2 Preliminaries and Theorem Statement 

Consider a two-dimensional, piecewise-C'^ continuous system of ordinary differential 
equations in with two independent parameters, and i]. For our analysis below 
we will need to assume A; > 8. In a neighborhood of a single switching manifold the 
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system may be written as 



X 

y 



f^^'>{x,y;fi,r]) 
g^^\x,y;fi,r]) 



Hix,y;n,r]) < 
H{x,y;fi,r]) > 



(2) 



where : M"* — M is a sufficiently smooth (at least C^) function. The switching 
manifold is the parameter dependent set, 5^^^ = {{x, yY \ H{x, y; /i, t]) = 0}. With- 
out loss of generality we may assume that H vanishes at (x, y; fj,, rf) = (0, 0; 0, 0). If 
in addition §f- + ^ 7^ 0, then locally iSo,o is a curve intersecting the origin. Via 
coordinate transformations in a similar manner as given in [H], we may assume to 
order that H is simply equal to x. The higher order terms in H do not affect 
our analysis below; thus in what follows, we will assume H is identically equal to x. 
The switching manifold is then simply the y-axis and we will refer to {f^^\ g^^^'Y as 
the left-half- system and {f^^\g^^^Y the right-half-system. 

We may assume that there is a discontinuous bifurcation at the origin when 
/i = ?7 = 0. Since the origin lies on the switching manifold and ([2]) is continuous, 
it is an equilibrium of both the left and right-half-systems. In this paper we are 
interested in the scenario that the equilibrium in the left-half-plane has complex- 
valued eigenvalues X± = p ± iu. (We make no assumptions about eigenvalues of 
the equilibrium solution in the right-half-plane.) Assume that when /j, = r] = 0, the 
eigenvalues are purely imaginary, i.e.. 



1/(0, 0) = 0, a;(0,0)>0 
Notice ±iu;(0, 0) are the eigenvalues of the matrix 



J 



Therefore, in particular 



dx 
dx 



det(J) 

Jl2 



dy 
dy 



(3) 



(4) 



(x,j/;/i,77)=(0,0;0,0) 



a;2(0,0) ^ 
dfiL) 



dy 



7^0. 



(5) 
(6) 



(0,0;0,0) 



By the implicit function theorem and ([5]) the left-half-system has an equilibrium 
{x*^^\jj,r]),y*^^\n,r])Y where x*^^^ and y*^^^ are C'' functions and a;*(-^)(0,0) = 
y*^^\0, 0) = 0. As is generically the case, we assume the distance of the equilibrium 
from the switching manifold varies linearly with some combination of parameters. 
Without loss of generality we may assume is a suitable choice. That is 



(0,0)^0. 



(7) 
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Again by the implicit function theorem, there is a function, 0, such that x**-^-* {(p{r]), r]) 
0. In other words when yu = ^(t]), the equihbrium hes on the switching manifold. 
After performing the nonlinear change of coordinates 



y ^ y - y*^^\(p{v),v) , 

we may factor fi out of the constant term in the system ([2]), i.eQ 



(8) 



/(^)(0,0;Ai,r/) " 




■ /(^)(0,0;/i,r7) " 




p{fx,r]) 


f7(^)(0,0;Ai,r/) _ 




_ (7(^)(0,0;/i,r7) _ 







fi + o{k) 



(9) 



where p and q are C . The left and right-half-systems may now be written as 



p(/i,77) 





X 







f'^'\x,y;fi,ri) 
^W(x,y;^,?7) 



+ o{k) 



(10) 



where Oj, b, Ci, d are each C''~^ and f^^\g^''^ are C'' functions nonlinear in x and y. 
(As with p and g, subscripts are not required for the coefficients b and d because the 
system is continuous.) 

We show now that it is possible to change coordinates to set p{fi, rj) = 0. Since 
f)(0,0) 7^ by (ini), using the implicit function theorem, there exists a unique C'' 
function, ip, such that f^^\0,tp{ix,r])] fj.,!]) = 0. After the coordinate change 



y^y- 



the system remains in the form (ITOl) with p(/i, rj) = 0. 

Finally we transform the system to the usual canonical companion matrix form. 
As is well-known, this may be accomplished when the system is observable in the 
control theory sense [Ul E]- Our system is observable by and the required 
transformation is 



y ^ -d{fM,r])x + b{fi,r])y , 
fi ^^ -b{fi,ri)q{fi,r])fi , 
7] t-^ aiifi,!]) + d{fi,ri) . 



(12) 



The first equation is nonsingular by ([6]). The second equation is nonsingular by ([7]) 
and the final equation is nonsingular if we make a final nondegeneracy assumption: 



3^(0,0)^0. 



(13) 



The transformed system has a normal form suited for further analysis and is given 
in the following theorem. 



^ We use 0{k) (and o(fc)) to denote terms that are order k (larger than order k) in all variables 
and parameters. When necessary to distinguish orders we are more specific, e.g. 0{\x,y\^). 
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Theorem 1 

Consider the two-dimensional, piecewise-C^ (k > 8), continuous system of dijferen 
tial equations 

{f^^\x,y;iJ,,r]) 
9^^\x,y;fi,r]) 
f^^^{x,y;^,r]) 
g^^\x,y;n,r]) 



x<0 
x>0 



(14) 



where 



f^'\x,y;^,ri) 
g^^{x,y;n,r)) 



such that 








+ 









1 " 




X 









y 


+ 



f^^{x,y;^,r]) 
f^{x,y;id,r)) 



+ o{k) , 
(15) 

(16) 



and f^'^Kg^^^ consist of terms that are nonlinear in x and y. Suppose, 
i) 5(^)(0,0) = for cu > 0, 



a) ao 7^ 0; where 
1 

ao -- 



\J XXX ' iJxxy ' J xyy ' '^yyy' 



16 

1 



-^qJ xy \J XX ' -lyy ■> 



16 ^(j^^^ 



yy ' ' 15^^2-'3;a; ^xx Jyy Hyy > 



(17) 



evaluated at {x, y; /i, i]) = (0, 0; 0, 0), and 

iii) Tr = r(^)(0,0) ^ 0. 

Then, near {x,y;fi,r]) = (0,0; 0,0), if Sr = 6^^\0,0) ^ when Sr/j, < there is a 
unique equilibrium in the right-half-plane given by functions 

Ai + 0(2) . 



x*'-^\lj,r]) 
2/*(«)(/x,r;) 



_ j_ 

Sr 

TR 

Sr 



(18) 



The equilibrium is repelling if tr, 6r > 0, attracting if —tr, 6r > and a saddle if 
5r < 0. When > 0, there is a unique equilibrium in the left-half-plane given by 
functions 



x*^^\n,r]) 




1 " 


_y*(^)(//,r/) _ 








l^ + 0{2) . 



Furthermore, there exist unique ^, C'^ ^, ^ functions hi,h2,hs 
respectively, with 



(19) 



h2{n) 
hsili) 

such that when fJ> > 0, 



1 



^2 ^-fxx^ fl'ij/O I (0,0;0,0)^ ) ' 



R 



(20) 
(21) 
(22) 
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i) the curve r] = hi{fi) corresponds to a locus of Andronov-Hopf bifurcations of 
{x*^^\ y*^^^)'^ that are supercritical if ao < and subcritical if ao > 0; this 
equilibrium is attracting if rj < and repelling if r] > hi{fi), 

a) the curve rj = /i2(/i) corresponds to a locus of grazing bifurcations of the asso- 
ciated Hopf cycle with the y-axis, 

Hi) if a^TR < 0, the Hopf cycle exists for values of rj between and h^i^fi), a 

periodic orbit of opposing stability exists for r] < hs^fi) if ao < and r] > h^^^) 
if ao > and the two orbits coincide at a locus of saddle-node bifurcations 
V = 

iv) if ooTji > 0, the Hopf cycle exists for r] > if ao < and r] < if 

ao > 0. 

Theorem [1] predicts essentially two different unfolding scenarios. These are illus- 
trated in Fig. [H 




Figure 1: Schematics showing unfoldings predicted by Theorem [T] when tr < 0, 
6ji > 0. In panel A, oq < 0, and in panel B, > 0. Double lines correspond to 
Hopf bifurcations or discontinuous bifurcations that create a periodic orbit. Dashed 
lines correspond to a collision of an equilibrium or periodic orbit with the switching 
manifold at which no bifurcation occurs. Solid lines correspond to saddle-node 
bifurcations of periodic orbits. Included are phase portraits showing local behavior. 
Stable orbits are colored blue; unstable orbits are colored red. Thin solid lines denote 
the switching manifold. 
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3 Example 



As an example consider the piecewise-C°°, continuous system 

2 1 

u = -a + —f3 + v + -u^ + , 
15 5 

5 13 1,^ .X ,1 1 , 

V = — a H — p u -\ { J — Ijv + \-u v\ . 

4 6^ 8 10^'^ ^ '8 10 ' 



(23) 



When a = (3 = 0, the origin is an equihbrium on a switching manifold. Its two 
one-sided limiting associated eigenvalues are =t:^i and — ^ ± -^i- This example 
exhibits the codimension-two scenario in which we are interested and satisfies the 
required non-degeneracy conditions. 

A bifurcation set for fl2^ is shown in Fig. [21 For small values of a and j3 the 
bifurcation set is a smooth distortion of Fig. [H panel B. However here there also 
exists an unstable periodic orbit. The predictions of Theorem [1] break down away 
from a = P = where this orbit collides with the stable orbit. The resulting 
locus of saddle-node bifurcations of periodic orbits intersects the saddle-node locus 
anticipated by Theorem [T] at a cusp bifurcation at {a,j3) ~ (0.019, —0.29). 



0.2- 



-0.2- 



-0.4 




-0.02 



0.06 



Figure 2: A bifurcation set of the example system ( |23|) . Sketches of representative 
phase portraits are included. The line styles and color scheme is the same as in 

Fig.m 



In order to compare the predictions of Theorem [T] and in particular the scaling 
laws (l20l) - (p2l) . we transform the system (l23l) to the form given in the theorem. The 
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switching manifold is the hne u = $v, therefore we let 



X = u 



(24) 



This particular example has been chosen because the transformation may be com- 
puted explicitly. We have 0(/3) = -§/?, y<^\(p{P) , (3) = and iA(a,/3) = 0. Com- 
bining the individual transformations ([HD, (ITT]) and ffT^ produces 



y = - — W-^)u + v, 

1 2 
/^=-T7t(^-/^)(«-77/3), 



V 



10' 2 



15 



In the transformed coordinates the system is 

1 












+ 


X 


-\ 




-/^ 




y . 









+ 



X 

y 



+ 



X 

y 



2 
+ 



V 
1 



where 



25x + 20y 
33 - 207] 

The values of the important constants are 



{ao,uj,TR,6R) 
(f{L) 

\Jxx ' ; I (0,0;0,0) 



25 1 11 
88'7f'~5'4 
250 80 



1089' 1089 



(25) 



X < 



X > 



(26) 



(27) 



(28) 



Fig. [3] shows comparisons of numerical computations of the curves hi, /12 and /is, 
with their predicted scalings. We find the numerical results are in full agreement with 
Theorem [H A calculation of error terms to quantitatively estimate the difference 
between the lowest order approximations and the true curves is beyond the scope of 
this paper. 



4 Proof Outline 

Here we present an outline for the proof of Theorem [T] that follows in §3 By 
assumption, when ^ > 0, there is an equilibrium, {x*^^\y*^^'^Y , in the left-half- 
plane. Close to this point there will exist an appropriate Poincare section, 11, in 
the left-half-plane, see Fig. |H The trajectory that begins from a point po on 11, 
spirals clockwise around and reintersects 11 at some point p^. We are 

interested in periodic orbits of the flow, thus when p^ = Po- 



8 




Figure 3: Numerical verification of Theorem [T] for the transformed example system 
fl26l) . Panel A shows hi{rj). Panel B is a plot of the square root of the difference 
between /i2(^) and hiij]). Panel C is a plot of the sixth root of the difference between 
h^{rj) and ^2(^)- The dashed lines in the three panels are the lowest order scaling 
predictions fl2U]) - fl22]) . respectively. 



It will be more convenient to use a second Poincare section, 11, that is a semi- 
infinite line intersecting the equilibrium and the origin. Artificially following the 
left-half-flow from po? we arrive at a point pi on 11. If pi were in the left-half-plane, 
by continuing from pi along the left- half-flow we would arrive at p^. However if pi 
is in the right-half-plane, we must first calculate a correction, following the left-half- 
fiow from a different point on 11, in order to arrive at p^. The correction results 
from the lack of smoothness at the switching manifold. The mapping between pi 
and p4 is known as the discontinuity map [HI [H] . 

To compute the discontinuity map, we follow the left-half-fiow back from pi until 
we arrive at a point, on the y-axis after a time Ti < 0. We then follow the right- 
half-fiow from p2 until the next intersection with the y-axis at ps after a time T2 > 0. 
Finally we follow the left-half-flow to p^ on 11 after a time T3 < 0. Let Vi, V2 and 
V3 denote maps relating to these three steps, respectively. The discontinuity map 
is then 

Pdm = P3 O ^2 O Pi . (29) 

As will be shown, Vdm has a |-type singularity. This is because the system is 
continuous on the switching manifold [HI [SI- When pi lies in the left-half-plane the 
discontinuity map is taken to be the identity map. 

A map from po to p^ may be derived by composing a map from 11 to 11 with Vdm 
and with a map from 11 to 11. However, a map from pi to Pq is simpler (where Pq is 
the point on H obtained by following the left- half-flow from p^) and equivalent. Let 
Pihf denote the map following the left-half-flow from H to itself and let T4 denote 
the corresponding transition time. Then pi is mapped to p^ via 

V = Vim o Vdm ■ (30) 

A fixed point of V corresponds to Pq = pi- Notice p^ = pi exactly when p^ = po. 

The most intensive computations in the following proof are in deriving the maps 
Vim, Vi, V2 and V3 described above. Once this is accomplished, V is obtained 
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Figure 4: Schematic illustrating a trajectory and Poincare sections relating to 
Theorem [1] when /i > 0. A solid [dotted] curve corresponds to following the left 
[right] -half- flow. The true trajectory is the thick curve. The Poincare maps con- 
structed in the proof are p2 = Pi(pi), Ps = ^2(^2), Pi = ^sbs), Pe = VmiPi)- 



by composition and it remains to derive its fixed points to find /i2 and obtain the 
function by locating where V has a saddle-node bifurcation. In order to derive 
expressions, not for /i2 and h^, but rather for /i2 — hi and — h2, we will intro- 
duce adjusted parameters rji and 772, that represent the deviation from hi and /i2 
respectively. 

5 Proof of Theorem [1] 



For ease of notation we expand the vector field in the left-half-plane to second order 
in X and y 



f^^\x, y; fj,,!]) = r]x + y + aix^ + a2xy + asy"^ + 0{\x, y\^) , 
g^^\x, y; jj, t]) = -^- 5^^'^x + bix^ + b2xy + b^y"^ + 0{\x, yf) , 

where coefficients vary with ^ and r]. 

Step 1: Compute equilibria, eigenvalues and hi. 

By the implicit function theorem, there exists a unique function. 



such that 



/(^)(x*(^)(^,r7),?/*(^)(/i,77);^,r/) 



(31) 



y*{L) 



0, for small rj. Via 
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substitution of a series expansion of a;*^^^ and y*^^^ into (I3T1) it is readily determined 
that 



X 



if^^V) = + 



hi_ ^ #^0,0 ) 



+ ^ — + 0(3) , 



where /i and r] subscripts denote derivatives. Near {x, y; fi, rf) = (0, 0; 0, 0), {x' 



(32) 



iL)^y*{L)y 



is the only equilibrium of the left-half-flow and by ( 1321) it is an equilibrium of the full 
flow (i.e. admissible), exactly when /i > 0. U 6^ 0, the corresponding equilibrium 
in the right-half-plane may be determined similarly. 

Let J^^^ denote the Jacobian of f|T^ when x < 0. For small /i and rj the matrix 



Jl 32 
k J4 



J(^)(x*(^)(^,r;),y*(^)(/x,r^);^,r^) 

has the complex conjugate eigenvalue pair 

A± = z/ ± i^ , 

where v and ^ are C^~^ functions of /i and rj. To first order we have 

By the implicit function theorem, there exists a unique C^~^ function hi : 
such that z^(/x, /ii(/u)) = for small fi. Furthermore 



(33) 



+ -r] + 0(2) . 



(34) 



hi{fi) 



2ai + 62 



thereby confirming (120|) of the theorem. Let 

r]i = ri- hi{fi) 

represent the deviation from the Hopf bifurcation curve, then 

K^,^i) = ^i(^ + o(i)). 



(35) 



(36) 



Step 2: Compute oq and introduce polar coordinates. 

Let vx^ denote the complex- valued eigenvector associated with X± for (155]) . In a 
standard manner, we construct a matrix P using the real and imaginary parts of 

1 



P 



-(ii -i4)^ -4i2i3 



which is well-defined (because j2 7^ for sufficiently small fi, rj) and non-singular. 
Let 



u 

V 



P 



-1 



X 

y 



x*^^\li,ri) 
y*'^^\lJ,v) 
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Then the left-half-system in {u, v) coordinates becomes 









V 







u 

V 



+ 0(2), 



where 

Letting 
then 



z = u + iv 



(37) 



i = A+^ + 0(2). (38) 

Following standard proofs of the Hopf bifurcation [TTl [TBI [19], there exists a two- 
variable polynomial comprised of only quadratic and cubic terms such that the 
near identity transformation 

u; = ^ + C(^,^), (39) 
removes all quadratic terms and all but one cubic term from the left-half-system 



where A G C and 
^e(A(/i,r/)) 



vj = A+w + Aw^w + 0(4) 



-t ^ i^fuUU ~l~ 9uUV ~l~ fuVV ~l~ Qvvv) ~l~ -t J. ifuvi^fuU ~l~ fvV^ 

lb 1d^ 



(40) 



9uv{,9uu ~l~ 9vv) fuuQuu ~l~ fvvQvv) 

Notice TZe{A{ii,r])) is a C^'^ function of and r/. Let 



(41) 



(n,i))=(0,0) 



ao = 7^e(A(0,0)) 



Using dH]) and since P(0, 0) 



1 

-u 



it follows that appears as given in the 



statement of the theorem, ffTTj) . 

To prove that generic Hopf bifurcations occur along ^ = hiij]) it remains to 
verify the non- degeneracy conditions of the Hopf bifurcation theorem (TTJ dSl (TU] . 

i) By construction, J'^-^)(a;**^^^(/i, /ii(yu)), y*^^^(/i, /ii(/i)); yU, /ii(/i)) has purely imag- 
inary eigenvalues, ±i^ = ±i(a; + 0{^)) ^ 0. 



dv 



(^,/ll(Ai)) 



0(/x) ^ 0. 



iii) By assumption, 7le{A{fj., /ii(/i))) = oq + 0(/i) 7^ 0. 
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Therefore if ao < 0, the curve i] = hi{fi) corresponds to supercritical Hopf bifurca- 
tions and stable periodic orbits exist for small r] > hi{fi). Conversely, if > 0, the 
curve T] = hi (/i) corresponds to subcritical Hopf bifurcations and unstable periodic 
orbits exist for small 77 < hi{fi). Consequently we have proven (i) of the theorem. 
We now introduce polar coordinates. Let 

w = re'^ . (42) 

In polar coordinates the left-half-system is 

r = ur + ne{A)r^ + ©(r'') , 

(43) 

e = C + Im{A)r^ + 0{r^) . 

Denote the components of the flow by R{rQ,9o,t] fj.,rii) and 6(ro, 6*0, t; /x, 771) 
respectively. Expressions for R and may be derived by expanding each as a series 
in To and computing coefficients by solving initial value problems. We obtain 

i?(ro,^o,t;/i,??i) = e-Vo + 7^e(A) + 0(4) , (44) 

e(ro,^o,t;/i,m) = eo + ^t + Im{A)^—rl + Oirl). (45) 

Step 3: Define the Poincare section, 11. 
By using fl3T|) it is straightforward to show 



lim 



y*{L) 



Thus the quotient 

g = £ 

is a well-defined function for small 77. For our analysis we only need the Taylor 
series of 5* to first order, which from fl32l) is readily found to be 

S{^l,r^) = ^^^^-n + 0{2). 

For /i > 0, ?/ = is the line intersecting (x*^^\ •')"'" and the origin. Let 

n = {(x,?/) I 2/ = > , (46) 

as in Fig. |H In the polar coordinates centered at the equilibrium, {r,9), U is de- 
scribed by a function 

9uir;^i,Vi) = + pir + O{r^) , (47) 

where 9 and pi are coefficients dependent on p and rji. The coefficient 9 describes 
the angle of 11 at the equilibrium in {u, v) coordinates; pi arises from the nonlinear 



13 



coordinate change, (l39l) . Since explicit forms for these coefficients will not be re- 
quired, we do not derive them. 



Step 4: Derive V\m and compute h2- 

We now wish to determine the left-half Poincare map. Vim- To do this we compute 
the trajectory of a point, {rQ,9-si{ro; ^jiji)) on 11, and find its next intersection, 
(n, Onii^i, f^, Vi))- The transition time is 

27r 

n{ro;^i,vi) = — + 0{l). (48) 

UJ 

Substituting fl36|) and psj) into fj44l) produces 

= (l + ^,(l + 0{\fi,vi\')))ro + ^rl + 0{4) . (49) 
Let s denote the distance from {x*^^\y*^^^) in {u,v), fl37|) . coordinates. That is 

II / 2 2 \ — 

s = |2;| = (tt + f j 2 . 



By (I39D and (02]) 

s(r) = r + (T2r'^ + O(r^) , (50) 

for some coefficient (T2, determined by (!39|) . By combining (H9|) and (l50l) we are 
able to obtain si as a function of sq, where sq is a point on 11 and si is the next 
point on 11 

si{so; /i, vi) =(l + Vi{- + 0(1//, ?7ir))) So + —Visl + + 0(4) . (51) 

Let {e,Se)'^ be a point on LI in (x,y) coordinates. Let 

e = e- . (52) 

Since e is a scalar multiple of s and when ^ = r]i = we have 5 = s, it follows that 
to third order the map between Eq and £1 is the same as fl5T]) . i.e. 



£i(£o; /X, ??i) = (1 + ?7i(- + 0(1//, Vi\'))) ^0 + — r/i£o + ^4 + 0(4) . (53) 

\ UJ / UJ UJ 

The system (IT^ has a periodic orbit that grazes the y-axis when 5o = £1 = —x*^^^ 7^ 
0. By substituting this into ( l53l) and dividing through by — x*'^'^) we obtain 



1 = l + vi{- + Oi\f,,v,\')) - !^r7ix*(^)(/x,r^i) + (/i, r^i) + 0(3) . (54) 

UJ UJ u 

Using fl32l) and the imphcit function theorem we find fl5^ is satisfied when r]i = /i2(/<) 
for a O'^ function 

2an 



14 



Let 



h2 

V2 



hi + h2 , 



(55) 



We have therefore derived fl2Tl) and proven (ii) of the theorem. We now write V\m 
in terms of e, /i and 772- Combining fl52|) . fl53l) and fl55|) produces 



'Pihf(£o;/^,??2) 



00 



where 



+ q3V2 + 0{\i^,V2n)eo + 0{e'o) , 
47rao 



(56) 
(57) 



and exphcit forms for q2 and ^3 will not be required subsequently. 



Step 5: Derive the discontinuity map, Pdm- 

In order to avoid singularities, here we assume fi > and introduce the spatial 
scaling 



X 


1 


X 









In the left-half-plane 



X 




. y . 





f'^^\x,y;iJ,,r]) 
g'^^\x,y;n,r]) 



1 
L 



f^\lJ,x,ny;^,vi) 
g^^\fix,fiy;fi,r]) 



Let 



X(^)(yo,t;^,r7) 

y(^)(yo,t;/^,r/) 



A,t + A^t" + ^3^3 + Oit") , 
yQ + Bit + B^t" + B^t^ + 0{t^) 



{51 



(59) 



(60) 
(61) 



denote the components of the left-half-flow for an initial condition (0, ijo), on the 
switching manifold. The parameter dependent coefficients are obtained by solving 
dSlD using ([SI]) 



A, 
A2 

A3 
Bi 

B2 
B3 



yo + a-iiiyl + 0{yl) , 

+ (^^ - 03/^)^0 + 0{yl) , 

-03^ - -r] + 0(yo) , 
3 6 

-1 + h3^^yl + O(yo') , 
(-^5^^) - &3/^)yo + O(yo') , 

lf,3^+l5W + 0(yo) . 

3 6 



(62) 
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We now derive Vi, V2 and V3 in scaled coordinates, (1581) . beginning with V3. The 
point p4 = Vzip^) (see Fig. Hj) hes on 11, therefore to compute the corresponding 
transition time, T3, we solve y(^) = SX^^") for t. The function 

= yo-t-Syot + ^St' + 0{3) , 

is and by the implicit function theorem there exists a unique function T3 : 
M'^ ^ R such that Gi{yo, T^O/o] /i, "'?) = 0. Furthermore 

T3{yo]fi,v) =yo-^Syl + 0{y^) . (63) 



Combining (l60l) and ( l63il yields 

^3(^0; ^, r?) = i^o' + + «3^)^o + 0{y^) . (64) 

Notice (see Fig. Hj) that Pf^ is the same as V3 except the y-component of the initial 
condition has opposite sign. That is, yo = Vi{io', /i, ?]), whenever £0 = 'Psi^Vo] A^i v)- 
By inverting flM|) we obtain 

Vi{io]fM,r]) = V2e^ - -{r] + a3n)io + 0{i^) . (65) 

In a similar manner as for V3 we are able to use a series expansion of the right-half- 
flow to determine V2- We obtain 

^2(^0; Ai, ^) = -yo - ^(t(^) + a3^)yo' + OiyI) . (66) 



The discontinuity map is Vdm = Vs 0V2 ^Vi. Composition of (l64l) - (!66l) produces 

A /o 3 

£4 = Vdmiii; /u, v) = Si + -^(^^""^ - + . (67) 

Step 6: Obtain the full Poincare map, V and compute /13. 

The full Poincare map is P = V\m o Vdm- In scaled coordinates fl56l) becomes 

^6 = r?2(^ + 0(|^, ^2r))+ (1 + + + q2^r]2 + 93^2 + 0(1^, r^a^) £4 + 0(£l) . 

(68) 

Composing (1671) and (168|) produces 



£6 = P(£i; /i, ??2) = f^o + + ^2el + 0(£f) , (69) 



where 



f^o(^,r/2) = ^r/2 + 0(2) , 

^i(m, ^2) = 1 + -??2 + giyw^ + g2Ai?72 + qzvl + 0(3) 
a; 



4^/2 

^^2(M,r72) = -^r^ + 0(l). 
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To remove fractional powers we introduce x = v^- The function 

G2(x; V2) = V{x^; /i, 772) - X^ , 

is C^~^, and by the imphcit function theorem, there exists a unique C^~^ function, 
JF, such that G2(x; A^j -^(xi /^)) = 0. Via a series expansion it is straightforward to 
obtain 



Hx; /i) = ( -^/i^ + 0(^^)1 + f -14^^ + o(^) I + 0{x') . (70) 



The fixed point ( iTOl) . has an associated multipher of one when the ^ function 

= (gi/i^ + 0(/i3)) + (2v^r« + + O(x') , 

is zero. By the imphcit function theorem, there exists a unique C^~^ function, h^, 
such that G^^h^^fj,); fj,) = 0. Furthermore, using (IFTIl . 

1 

Notice this fixed point is vahd when e\ = x = ^3(/^) ^ which is true when ao and 
T/j have opposite signs. Finally let 

Then is the C^'"^ function (1221) and we have verified (iii) and (iv) of the theorem. 

□ 



6 Conclusions 

We have presented an unfolding of the codimension-two simultaneous occurrence of a 
discontinuous bifurcation and an Andronov-Hopf bifurcation in a general, piecewise- 
smooth, continuous system. We have found a locus of Hopf bifurcations that em- 
anates from the codimension-two point, fl20|) . Tangent to this is a locus of grazing 
bifurcations of the Hopf cycle with the switching manifold, ( !2T|) . Heuristically, the 
curves are tangent because, with respect to a linear change in parameter values, 
the distance between the equilibrium solution and the switching manifold increases 
linearly, whereas the amplitude of the Hopf cycle grows as the square root of the 
magnitude of the parameter change. 

A periodic orbit is created at the discontinuous bifurcation on one side of the 
codimension-two point. When the stability of this orbit opposes that of the Hopf 
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cycle, the two orbits collide and annihilate in a saddle-node bifurcation on a curve 
that deviates only to order six from the grazing bifurcation (J22i) . The mechanism 
behind this sixth order scaling can be explained with a simple calculation. Omitting 
higher order terms, the Poincare map, "P, fl69l) . is essentially 

e' = 772 + H(/i)£ + 7£i , (71) 

where is the Floquet multiplier of the Hopf cycle and 7 is a constant. The 
grazing bifurcation occurs when 772 = 0. It is easily determined, as in [8J, a saddle- 
node bifurcation of the fixed point of ( ITTI) occurs when 



The Floquet multiplier is unity at the Hopf bifurcation (t/i = 0) and by assumption 
varies linearly with respect to r^i, i.e., S ^ 1 + Ar^i, (A 7^ 0). The grazing bifurcation 
occurs when r^i = 0(/i^), thus S 1 + A/x^, (A 7^ 0). Hence 

For simplicity, throughout this paper we assumed the switching manifold was 
infinitely different iable. Our analysis is unchanged if the switching manifold is only 
C^. However, if the switching manifold were and not C^, we would be unable to 
determine the same expression for the map P3, fl64l) . 

Recently it has been found that an eight-dimensional model of yeast growth |T3] 
exhibits codimension-two discontinuous bifurcations such as the scenario described 
here. Other observed codimension-two situations that remain to be rigorously un- 
folded include the simultaneous occurrence of a saddle-node and discontinuous bi- 
furcation, and a discontinuous Hopf bifurcation P, [TU] of indeterminable criticality. 
We hope to report on these in a future paper, see also [20] . 

We would also like to extend the results of this paper to higher dimensional 
systems, like the yeast model. It seems plausible that bifurcation sets for higher- 
dimensional systems will exhibit scalings of the same orders, but we do not, as of 
yet, have a formal justification of this. 
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